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In this work, the ratchet dynamics of Brownian particles driven by an external sinusoidal 
(harmonic ) force is investigated. The gating ratchet effect is observed when another harmonic 
is used to modulate the spatially symmetric potential in which the particles move. For small 
amplitudes of the harmonics, it is shown that the current (average velocity) of particles exhibits 
a sinusoidal shape as a function of a precise combination of the phases of both harmonics. By 
increasing the amplitudes of the harmonics beyond the small-limit regime, departures from 
the sinusoidal behavior are observed and current reversals can also be induced. These current 
reversals persist even for the overdamped dynamics of the particles. 

PACS numbers: 05.40. -a, 05.45. -a, 05.60.-k 


1 Introduction 


The transport of particles or solitons under zero-average forces (i.e., ratchet transport) has been extensively 
investigated in the last two decades mm mm in]. This phenomenon has been predicted and explained in 
different fields of physics, ranging from nano-devices to molecular motors Hamuli]- Moreover, it has also been 
observed in experiments and simulations with nonlinear systems, where spatio-temporal symmetries have been 
properly broken In particular, the ratchet models were used: to elucidate the working 

principles of molecular motors; to design molecular motors H31; and to explain the unidirectional motion of 
fluxons in Josephson junctions [23JI2n]> the transport of cold atoms in optical lattices [M], and the vortices in 
superconductors Earn. 

The ratchet transport is described by means of the current (average velocity) El CCS HU HI] ; 


v = lim 

£—>■00 


(x(t)) - Xp 
t-t 0 


(1) 


where x(t) is the position of particles, or the center of mass of solitons at time t , (•) represents an ensemble 
average over all trajectories satisfying the same initial condition, and x{to) = Xq. 

Two possible underlying mechanisms of rocking ratchets are harmonic mixing and gating. The current of 
particles (atoms or solitons) in harmonic mixing is generally induced by an additive bi-harmonic, T periodic, 
driving force /(f) = /Rt) + / 2 (t), with 


/i(t) = ei cos(giwt), / 2 (t) = e 2 cos(g 2 wt + />), 


( 2 ) 


where e\ and e 2 are the amplitudes of the harmonics, <f> is the relative phase between the two harmonics, 
(gi,g 2 ) £ N 2 , gcd(gi,g 2 ) = 1 and T = 2-k/uj. On the other hand, in gating ratchets , particles experience a 
symmetric potential with the amplitude modulated by means of fi(t). A time-symmetric harmonic force / 2 (t) 
is also applied. 
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The time-shift invariance of the current, 


(3) 


v[fi(t + r),f 2 (t + T)] = v[f 1 (t),f 2 (t)}, 


Vt, together with the symmetry 




(4) 


or 

v[fi(t),-f 2 {t)] = -v{fi{t),f 2 (t)\, (5) 

fix the necessary conditions on q\ and q 2 in Eq. ([2]) to obtain the ratchet effect in harmonic mixing and gating. 
Symmetry Q holds for rocking ratchets induced by an additive bi-harmonic force, whereas ([5| characterizes the 
gating average velocity. When q\ +q 2 is an odd integer number, the bi-harmonic force /(f) breaks the time-shift 
symmetry /(f) = —f(t + T/ 2) and a current appears. In a gating ratchet, if q\ is an odd integer number, /i(f) 
preserves the time-shift symmetry, where /i(f) = —/i(f + T/2). Nevertheless, the gating effect appears due to 
a synchronization of the oscillations of the potential barrier caused by a single harmonic /i(f) with the motion 
produced by the additive harmonic force, / 2 (f)- There is no constraint on q 2 in gating, and therefore a current 
can be obtained even for q\= q 2 = 1 mm- 

Moreover, Eqs. <M together with the functional representation of the ratchet velocity determine the 
dependence of the current on the amplitudes and relative phase of the harmonics mm- For instance, for the 
small-amplitude limit of the bi-harmonic force /(f) = /i(f) + / 2 (f) with ([2]), the current reads 

v[f{t)] = A 0 efef cos(<?i</ + 0 O ), (6) 


where q\ + q 2 is an odd integer number. Otherwise the current v vanishes mm- The constants Aq and 9q are 
determined by the other parameters of the system (potential, dissipation, etc). Equation ([6]) clearly shows the 
harmonic mixing since the parameters of the fist harmonic always appear in combination with the parameters 
of the second harmonic. Interestingly, for a gating ratchet, it is deduced (for a small-amplitude limit) that v 
again is ruled by Eq. however only q\ should be an odd integer number, whereas q 2 can be either an odd 
or even integer number. This formula predicts a sinusoidal dependence of v versus the phase (f>. This implies, 
for example, that current reversals can be induced by solely changing the relative phase between /1 and f 2 . 
Furthermore, in )5], for a non-snrall amplitude limit, two interesting effects have been theoretically predicted: a 
deviation from the sinusoidal shape of v as a function of the phase; and the dependence of 0 O and A 0 on the 
amplitudes of the forces. This latter fact leads to an unexpected phenomenon related with the appearance of 
current reversals by changing the amplitudes of the harmonics. This explains the experiments in optical lattices 
driven by a bi-harmonic force reported in [3j, and in a shaken liquid drop driven by two independent harmonics 

imi- 

iii this work, we focus on the ratchet dynamics of Brownian particles lying in a symmetric potential, modulated 
by a harmonic function. The particles are driven by an external sinusoidal ( harmonic ) force. We show that there 
is a deviation from the sinusoidal behavior of v as a function of the relative phase between the two harmonics 
in the non-small amplitude limit. Moreover, the current reversals by means of increasing the amplitudes of the 
harmonics are shown. 

The paper is organized as follows: In the next Section, the symmetry properties of the Langevin equation 
and its relation with the functional representation of the current predicted in [5] are described. In Section III, 
the analytical predictions of the previous section are verified by means of simulations. In addition to a class of 
current reversals, determined by dissipation-induced symmetry breaking [3], we show that the current reversals 
persist even for the overdamped dynamics of our model. To conclude the paper, in the last Section, the results 
of Sections II-III are discussed, thereby making the connection with the experiments and summarizing our main 
findings. 
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2 Gating ratchet model 

In our theoretical analysis, the dynamics of particles in the spatially symmetric potential is determined by the 
Langevin equation 

mx = -ax - U\x)[l + /i (t)] + f 2 (t) + V2D^(t), (7) 

where m is the mass; U{x) = UqCOs(x) is a periodic symmetric potential, modulated by the harmonic fi{t) 
given by ([2]); a the friction coefficient; £(i) a Gaussian white noise, (£(f)) = 0, = 8{t — t'). Generally, 

noise smooths the dependence of the current on the parameters of the harmonics [18] . In some cases, as we show 
below, adding noise promotes transport. The additive force f 2 (t) are given by ([2]). All these magnitudes and 
parameters are in dimensionless form. 

The current defined by Eq. 0 is time-shift invariant, i.e. it fulfils the symmetry 0 due to the dissipation. 
Therefore, if v is a smooth functional such that its functional Taylor series exists, then Theorem 1 of [5] assures 
that 

OO 

v = YXcrWC&u £ 2) COS (kq 1 (f> + 0fc(ei,e 2 )), (8) 

k=0 

with 0 o(ei,e 2 ) = 0, and functions Cfc(ei,£ 2 ) and the phase lags 0fc(e 1 , 62 ) are even in each e 7 , j = 1,2. Notice 
that the symmetry ([ 5 ]) holds since exchanging f 2 with —f 2 is equivalent to replacing x(t) with —x(t) in (J7|. The 
statistical properties of the Gaussian white noise are the same under the inversion of £(f) to — £(i). Therefore, 
all Ck with even fc are zero. With this restriction, the first two terms in ([8]), for qi = q 2 = 1, read 

v = vicos(<j) + 9i) + v 2 cos(3^ + 0 2 ) + Ei 0 (ei,e 2 ), (9) 

where v\ = eie 2 Ci(ei,e 2 ), v 2 = (ei£ 2 ) 3 C l 2 (ei, £ 2 ), C 1 and 9\ are polynomials up to order 6 in e 7 , and C 2 and 9 2 
are linear in e\ and e^. 

For (71 = 1 and q 2 = 2, v is given by 


v = vi cos {<{> + 9\) + v 2 cos(3 (j) + 9 2 ) + Si 5 (ei,e 2 ), 


( 10 ) 


where E 15 contains terms of order 15 or higher in each ey, v\ = efe 2 Ci(ei : e 2 ), v 2 = (ef e 2 ) 3 C 2 (ei, e 2 ), C\ and 9\ 
are even polynomials in and £2 up to order 10, and C 2 and 9 2l are even polynomials in £i and £2 up to order 
4. In both cases, (q± = q 2 = 1 or q-\ = 1 and q 2 = 2) we have identified 3 main regimes which depend on the 
amplitudes of the harmonics, namely: 


1. Small-amplitude regime. Only the first term in ([9| and (10) dominates and C\ and 9\ do not depend on 
the amplitudes. By fixing all the parameters of the system, v is a sinusoidal function on <j). 


2. Intermediate amplitude regime. The second term in 0 and ( |Toj ) can be neglected. However, in contrast 
to the previous case, C\ and 9i do depend on ei and £ 2 - Therefore, the current reversals can be achieved 
by modifying the amplitudes. The sinusoidal behavior of v persists. 


3. Large amplitude regime. The effect of the second term in 0 and ( fl()| ) is observed, and therefore v is no 
longer a sinusoidal function. 


The previous analysis remains valid for the overdamped dynamics. To describe the overdamped system we 
set m —> 0 in Eq. (|7|) : 

ax = -U'(x)[l + fi(t)] + f 2 (t) + V 2 D£(t). (11) 
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Moreover, time-reversal now implies that by changing fi(t) —>• f 2 (t) —> f 2 {—t) and x(t) —► tt — x(—t), 

the Eq. © remains invariant and 

( 12 ) 


This symmetry fixes all the phase lags in Eq. ([8]) to zero. Therefore, all the phase lags in Eqs. © and flip] ) are 
also zero. Nevertheless, current reversals may still be observed by changing the amplitudes of the forces. For 
instance, in the intermediate regime, a variation in the parameters around the values for which vi(ei, e 2 ) = 0 in 
Eqs. |9]) and (101, could make v change its sign. 

In the following section, all these findings are verified by means of simulations of Eqs. ([7]) and ©■ 


3 Simulations of the Langevin equation 

Simulations of the stochastic differential Eqs. 0 and © have been performed using the Heun method and the 
2nd-order weak predictor-corrector method P3]. The final time of integration is 2000, the time step is either 0.1 
or 0.01, and results are averaged over 10000 realizations unless specified otherwise in the figure caption. 





Figure 1: v vs (f> from simulations of ([7|. Filled circles with error bars: current in steady state computed from 
the Eq. (13). Solid line represents the fitting curve of the circles, v = —0.00170 cos(cj) + 2.195). Parameters: 
to = 1, a = 1, Uq = 5, ei = e 2 = 0.5, to = 1, q\ = 1, q 2 = 2 and D = 1. 


The current v is computed by means of 


v = 


x (tf) x{t r ) \ 
tf — t r / 


(13) 


where t r and tf are the final time of integration and the transient time, respectively (see Fig. [I]). 

The sinusoidal behavior of v((j>) is characteristic of the small and intermediate amplitude regimes. In Fig. 
HJ a sinusoidal behavior of v is observed as a function of the phase. Close to </> « 2.5 and </> ~ 5.7, the velocity 
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Figure 2: v vs <f> from simulations of (|7|) shows non-sinusoidal behavior (filled circles with error bars). Solid line 
represents the fitting curve v = — 0.4168cos(</> + 0.7713) — 0.0686cos(3<^> — 0.0909). Parameters: to = 1, a = 1, 
Uq = 5, ei = e 2 = 2, oj = 1, qi = 1, q 2 = 2, and D = 1. Final time of integration 1000. 


changes its sign and current reversals can appear by varying the phase and other parameters of the system that 
have an influence on the phase lag. 

By further increasing the amplitudes, the average velocity deviates from purely sinusoidal behavior and 
sinusoids of higher frequencies appear in its expansion. Indeed, in Fig. [2j the results from simulations of Eq. (J7| 
can be fitted perfectly with two harmonics. In Figs. (|T|) and Q, we notice that on replacing (j> with </> + 7r (this 
is equivalent to replacing / 2 with — / 2 ), v changes its sign. This means that the symmetry ([5]) is fulfilled. 

By fixing all the parameters, except ei and e 2 which vary according to e± = Ae, e 2 = A( 1 — e), we verify that 
the dependence of v on e is different from the expected v ~ e 2 (l — e), which is valid for the small-amplitude 
regime (see Fig. [3|. 

In order to observe a current reversal via an amplitude change, first we fix all the parameters of Eq. ([T]) as in 
Fig- [j] except the amplitudes of the harmonics, which we have increased up to ei = e 2 = 1. The amplitudes are 
now sufficiently large for the phase lags 9k to be no longer constant and for them to depend on the amplitudes 
ei,e 2 as in Eq. ([8]). We set a relative phase (j) « 2.8 which corresponds to an almost vanishing current for 
ei =€2 = 1 (not shown in the figures). A clear current reversal appears by modifying only the amplitudes 
around these values following e\ = 2e, e 2 = 2(1 — e) with e £ [0,1], as shown in Fig. [4j The inversion of the 
current occurs around e = 0.5, which corresponds to values = e 2 = 1 and a vanishing v, as expected. 

3.1 Overdamped dynamics of Brownian particle 

Interestingly, the maximum current shown in Fig. [5] for the overdamped particle is greater than the maximum 
current reached when the inertial term remains in the Langevin equation, see Fig. [I] Notice that the parameters 
in both figures are the same, except the inertial term which is omitted in the simulations reported in Fig. [5] This 
effect resembles the enhancement of the movement due to the dissipation studied in [221 EB hr the relativistic 
particle driven by a bi-harmonic force. 

This striking phenomenon vanishes when the amplitudes are increased (the maxima of the currents shown in 
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Figure 3: v vs e from simulations of 0 - Filled circles: current from the Eq. ( p~3| ) . Dashed line: fitting curve 
v = — 0.0207e 2 (l — e), predicted for small-amplitude limit. Solid line represents the fitting curve of the simulation 
points predicted in the intermediate amplitude regime, v = Ylk=o a kC k+2 i with a o = —0.002, aq = —0.076, 
a 2 = 0.367, 0,3 = —0.747, 04 = 0.664, 05 = —0.205. Parameters: m = 1, a = 1, Uq = 5, ei = Ae, €2 = A( 1 — e), 
A = 1, (j> = 4.09, u) = 1, qi = 1, 92 = 2, and D = 1. 



Figure 4: v vs e from simulations of ([T]) shows that the direction of the current changes at approximately e = 0.5. 
Parameters: m = 1, a = 1, Uq = 5, ei = 2e, e 2 = 2(1 — e), oj = 1, </> = 2.8, qi = 1, q 2 = 2, and D = 1. Final time 
of integration 1000. Dotted line represents zero velocity. 


Figs. [ 2 ] and [6] are almost the same). On increasing the amplitudes, a small deviation from the sinusoidal behavior 
of v as a function of the phase (j) is also observed in the overdamped system, see Fig. [6] 
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<l> 


Figure 5: v vs (f> from simulations of the overdamped system (11) (filled circles with error bars), 
represents the fitting curve v = 0.0033 cos(</>). Parameters: a = 1, Uo = 5, ei = £2 = 0.5, uj = 1, q\ = 
and D = 1 . 


Solid line 
1 ) 92 = 2 , 



0 1 2 3 4 5 6 

<t> 


Figure 6: v vs (f> from simulations of the overdamped system ( |11| ), filled circles with error bars. Dashed and solid 
lines are the fitting curves v = —0.345cos(0) and v = —0.344cos(^) — 0.042 cos(3</>), respectively. Parameters: 
a = 1, Uq = 5, ei = £2 = 2, 00 = 1, q\ = 1, q 2 = 2, and D = 1. 


In the overdamped dynamics, the phase lags are fixed to zero and the search for current reversals associated 
to changes in the amplitudes of the harmonics becomes a more difficult task (see Fig. [ 5 ]). In order to observe a 
current reversal by changing the amplitudes of the forces, we must proceed in a different fashion. Results from 


7 





simulations shown in Figs. [5] and [b] reveal that, by changing the amplitudes £i and £2 from 0.5 to 2 when <fi « 3, 
the direction of motion can be inverted. Therefore, by setting the phase, for instance at <j> = 2.8, and varying 
the amplitudes in the form of t\ = 2e, £2 = £1 with e £ [0,1], an inversion of the current is expected for a value 
of £i between 0.5 and 2. These results are shown in Fig. [7] Finally, Fig. [8] shows that an inversion is also 
observed when the amplitudes are modified while keeping the total amplitude £1 + €2 constant. Moreover, Fig. 
[8] shows that v = 0 when e = 0 (no modulation of the potential) or e — 1 (no additive force). 

In Figs. [T] and [8j the inversion of the current occurs at £i ~ 1.5. Indeed, by fixing all the parameters and 
changing and £ 2 , the contour plot (left panel in Fig. [ 9 ]) shows that the current vanishes when £i ss 1.5. 
However, for other sets of parameters, for instance taking Uq = 2.5 (see right-hand panel of Fig. [9]), the reversal 
current appears for different values of e\. 
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Figure 7: o vs e 
Amplitudes of f\ 
U 0 = 5, u = 1, qi 


from simulations of (11) shows that the direction of the current changes around e ~ 0.75. 
and /2 are varied simultaneously as e± = 2e, £2 = 2e. The rest of parameters are: a = 1, 
= 1, q 2 = 2, (j> = 2.8, and D = 1. 20000 realizations. 


4 Summary 

In this work, we study the dynamics of particles, driven by a harmonic force and subjected to white noise, when 
they are in a spatially symmetric potential that is modulated by a harmonic function. Both the applied force 
and the modulation of the potential are time symmetric and the current v fulfils the symmetry <©• Dissipation 
is also included in the description; therefore the current is time-shift invariant and the theory developed in [17 
can be applied. 

We show that this theory predicts three different regimes for our system which depend on the amplitudes of 
the two harmonics, namely: i) A small-amplitude regime where the current, v ~ e\ 2 cos(0 + $i), is a sinusoidal 
function with a phase lag, 9 \, independent of the amplitudes. This regime has been predicted by the collective 
coordinate theory and confirmed by simulations in the framework of soliton ratchets (see [57] and references 
therein), ii) The intermediate amplitude regime, where v is still a sinusoidal function, although 6 \ is no longer a 
constant. This means that the sinusoidal behavior alone cannot guarantee that the amplitudes of the harmonics 
are small. Therefore, in addition to experiments on optical lattices reported in [8], in order to determine the 
regime where the system lies, it is necessary to investigate the dependence of v on the amplitudes £i and £ 2 . Once 









Figure 8: v vs e from simulations of © shows that the direction of the current also changes around e ~ 0.75 
(filled circles) when the amplitudes are varied with constant total amplitude 2 as ei = 2e, e 2 = 2(1 — e). The 
rest of parameters are a = 1, Uq = 5, oj = 1, q\ = 1, qi = 2, (j) = 2.8, and D = 1. 20000 realizations. 




Figure 9: Contour plots v as functions of and e 2 from simulations of ©• The thick blue line marks the 
inversion of the current. Left panel Uq = 5. Right panel: Uq = 2.5. The rest of parameters are a = 1, u> = 1, 
<7i = 1, <72 = 2, (j) = 2.8, and £1 = 1. 


the intermediate regime is reached, current reversals via an amplitude change is expected. This phenomenon is 
confirmed by simulations of the underdamped Langevin Eq. |7]) . It is worthy of note that current reversals have 
been found to be present in the overdamped limit, where the current satisfies the time-reversal symmetry (12). 
iii) Large-amplitude regime, where we show that the non-sinusoidal behavior of the current, predicted by the 
theory, is due to the increasing strength of the two harmonics. 

Apart from the results presented in Figs.1-8, we have also performed simulations for all the set of parameters 
of Figs. l|8 but we fixed the strength of the noise D = 0. In all cases, the computed current is zero (of order 
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of 10~ 9 or less). Therefore, for the set of parameters studied here the noise together with the action of the 
harmonics generate the transport. 

Finally, it is pointed out that, according to the theory developed in [5], the main phenomena studied here using 
a specific model, can appear in other physical systems that satisfy the same symmetries, including experimental 
realizations in Josephson junctions 01251 eg and optical lattices [ 23 ], in which a number of the above results 
have been reported. Other results, however, require verification through experiments. 
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